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Abstract —We present an in-situ antenna characterization 
method and results for a “low-frequency” radio astronomy 
engineering prototype array, characterized over the 75-300 MHz 
frequency range. The presence of multiple cosmic radio sources, 
particularly the dominant Galactic noise, makes in-situ char¬ 
acterization at these frequencies challenging; however, it will 
be shown that high quality measurement is possible via radio 
interferometry techniques. This method is well-known in the 
radio astronomy community but seems less so in antenna mea¬ 
surement and wireless communications communities, although 
the measurement challenges involving multiple undesired sources 
in the antenna field-of-view bear some similarities. We discuss this 
approach and our results with the expectation that this principle 
may find greater application in related fields. 

I. Introduction 

As with any antenna engineering project, characterization of 
low-frequency radio astronomy arrays (referred to as “aperture 
arrays” to distinguish them from dish antennas) is required to 
validate that the design meets requirements and that modeling 
tools involved in the process are reliable. Low-frequency 
aperture arrays (LFAA) are particularly challenging to measure 
for a combination of reasons: they are difficult to point 
mechanically, are more closely coupled to the environment 
(see array placement in Fig. [T]) and their physical area is 
large for a given directivity (scales by A 2 ). As a consequence 
of the latter, larger distances (scales by the square of the 
array diameter) are required to satisfy the far-field condition. 
Without specialized facilities, measurements of LFAAs in an 
anechoic chamber or outdoor range are not practical Q The 
capability of testing larger-diameter LFAAs will be important 
in the pre-construction work of the Square Kilometre Array 
(SKA) radio telescope [2). 

In-situ LFAA measurements offers a number of obvious 
advantages: no bespoke facility is required and the array is 
situated in the intended environment. Fig. [I] illustrates the 
AUT: an array of 16 dual-polarized log-periodic “SKALA” 
antennas j3) deployed at the Murchison Radio-astronomy 

Unmanned aerial vehicles enable outdoor measurements 111, however the 
airborne height required to meet the far-field condition limits the size of the 
antenna or array under test (AUT). 


Observatory (MRO) in the Mid West region of Western Aus¬ 
tralia ©• Unlike measurements in an anechoic chamber where 
the operator has full control of the RF sources, astronomical 
sources are beyond our control. This is especially true in 
low-frequency radio astronomy, where the spatially extended 
component of the emission from the plane of the Milky 
Way Galaxy is the dominant source of noise 0. Meaning¬ 
ful LFAA measurement, therefore, requires exclusion of the 
diffuse Galactic emission while maintaining responsiveness to 
previously well-characterized, point-like extragalactic sources; 
this may be accomplished by cross-correlating the output 
voltages of pairs of AUTs which are spaced with sufficient 
distance such that the Galactic noise is uncorrelated. This 
concept will be discussed in Sec. [IT] Measurement results are 
and conclusions are presented in Secs. [TTT| and [Tv| respectively. 



Fig. 1. An LFAA prototype array of 16 SKALA antennas. This system is 
referred to as Aperture Array Verification System 0.5 (AAVS 0.5), and is an 
initiative of the SKA Array Design and Construction Consortium (AADCC). 
It was constructed by the International Centre for Radio Astronomy Research 
(ICRAR), University of Cambridge and the Netherlands Institute for Radio 
Astronomy (ASTRON). The SKALAs are placed in a pseudo-random config¬ 
uration within an area 8 m in diameter. 









II. Background Theory 
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A. Measurand 

The sensitivity of a radio telescope is quantified by the ratio 
of its effective aperture area (A in m 2 ) to the system noise 
temperature (T sys in K): 

a/t = A- (d 

1 sys 

This figure-of-merit is a key term in the expression for the 
telescope output signal-to-noise (S/N) ratio; many readers 
will recognize that it is similar and is convertible to G/T 
which is common in telecommunications. The antenna beam 
pattern may be obtained by measuring A/T along the apparent 
trajectory of a suitable astronomical source. This assumes that 
the sky noise does not change T sys appreciably during the 
observation; in practice, this is a reasonable approximation. 

B. Measurement Method 

This subsection presents a simplified version of how A/T is 
obtained in radio interferometry. Interested readers are referred 
to standard textbooks such as © and 171 for more complete 
discussions. Let us assume an interferometer involving 3 
single-polarized antennas (or antenna arrays as the case may 
be) observing a point source with flux S (in W/m 2 /Hz or 
Jy) in its phase center (in this context “phase center” refers 
to that position in the sky for which the interferometer is 
intended to point to by applying relative phase shifts between 
its elements). The measured cross-correlation products are: 

V12 = did^S 
V13 = did^S 

V 23 = a*atS , (2) 

where a is referred to as complex “gain” (voltage quantity) 
containing a combination of both antenna and receiver elec¬ 
tronic gains. Assuming S is known, the complex “gains” may 
be solved from these 3 equations. This is typically done by 
solving the amplitudes and phases separately ©• 

To obtain antenna gains, amplitude-only solutions are suf¬ 
ficient. Note that this method is analogous to the well- 
known “three-antenna” measurement technique With N 
elements, N(N — l)/2 pair-wise combinations (“baselines”) 
are available, leading to an overdetermined system solvable 
with a least-squares method. As shown in Fig. [2] the AAVS 0.5 
is co-located with the Murchison Widefield Array (MWA) 
radio telescope 0. PI to take advantage of this multi¬ 
baseline arrangement (with N = 128) to minimize calibration 
uncertainties. Each of the MWA’s 128 “tiles” is a 16-antenna 
array of bow-tie dipoles. 

By selecting interferometry baselines larger than a certain 
minimum, the dominant spatially-extended Galactic noise is 
no longer correlated and “bright” point-like sources are de¬ 
tectable. This may be illustrated by examining the spatial 
coherence function due to a distributed source (Ch. 1 in [6j): 

V 12 = [ [ I(l,m)e- j2 ^ Ul2l+Vl2m) dldm, (3) 
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Fig. 2. A map showing the AAVS 0.5 location (blue circle) and MWA 
telescope tiles (black cross) in Easting/Northing for Map Grid Australia UTM 
zone 50. 


where V± 2 is the complex visibility on the baseline between 
antennas 1 and 2, /(/, m) is the sky brightness distribution, u, v 
are the baseline vector components normalized to wavelength 
and /, m are the direction cosines with respect to v of a unit 
vector s = lu+mv + V 1 — l 2 — m 2 w where w is a unit vector 
normal to the uv plane (i.e., points to the “phase center”). 

Simplifying further to the 1-D case, Q becomes 

v 12 = J I(l)e- j2n(ui2l Ul, (4) 

where l = sin#' and O' is the angle between w and s. We 
recognize V\ 2 and 1(1) as a Fourier transform pair. Hence, 
structures with large l extent are mostly visible for small ui 2 . 
From MWA observing experience, with its small-diameter 
(3.3 m) tiles, correlated Galactic noise is measurable on base¬ 
lines of length < 30A and should therefore be excluded]^] 

It can be shown (see Ch. 9 of (6)) that after calibration of 
the measured cross-correlation products, the standard deviation 
in flux density for each baseline is inversely proportional to 
the square root of the products of the A/T for the AUTs in 
question: 


A&. 


k 2 (A/T)i(A/T)jBt a 


( 5 ) 


where k is Boltzmann’s constant, B is the bandwidth and 
£ acc is the accumulation time. Thus, by measuring AS for all 
baselines for which the Galactic noise is uncorrelated, A/T 
for each AUT may be solved. 


III. Characterization Results 

To measure A/T through cross-correlation, we point the 
interferometer towards Hydra A, a compact calibrator source 
that is sufficiently bright for a good S/N calibration solution 
to be obtained from a 2 minute snapshot observation (Fig. [3]). 
The AAVS 0.5 and MWA tile pointing is achieved with 5-bit, 
time-delay analog RF beamformers. The data is channelised 
prior to correlation into channels of width B -40 kHz, and 
averaged post-correlation to t acc =4 s. For each channel, A Sij 
is calculated for all samples in the snapshot period. 


2 Some readers will recognize a similar situation in a highly multipath 
environment: the signal envelope is uncorrelated at locations separated in 
distance by >~ 1A GU 
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Fig. 3. X-polarization cross-correlation image of Hydra A at 119 MHz, using 
only baselines between MWA tiles and the AAVS 0.5 array. The red ellipse 
(bottom-left corner) shows the size of the cross-correlation beam. 


A. Sensitivity 

Fig. [4] shows an A j /T measurement for a snapshot, made on 
22 May 2014, where the channels comprising the 30.72 MHz 
bandwidth have been grouped at 6 points across the available 
75-300 MHz frequency range. The top panel shows the X (E- 
W) antenna polarization in cyan dots, the lower panel shows 
the Y (N-S) polarization in blue dots. In this observation, 
both the AAVS 0.5 array and the MWA tiles were pointed at 
Hydra A, at azimuth (clockwise from North) Az=0°, elevation 
El=75.4°. The scattering of data points at ^270 MHz is due to 
the persistent satellite-based RFI corrupting the measurement 
in these channels. 

Fig. [4] also shows simulated results of this measurement as 
dashed curves; these follow Q and include measured receiver 
temperature, and the sky model ED and array gain pattern 
at the time of observation. We used full-wave electromagnetic 
(FEKC^ simulation to determine the array gain pattern for 
each polarization, and with a ground of either 2% or 10% 
moisture to represent a reasonable range of the likely soil 
moisture level experienced at the MRO | [I3) . 

The measured and simulated results show good agreement 
at all frequency points. With the exception of the corrupted 
data at ^270 MHz, the tight clustering of the results at 
each frequency point gives us confidence in the measurement 
process. 

B. Beam pattern 

Sensitivity can also be measured at Az, El angles away 
from the beam pointing direction enabling investigation of 
beam pattern. This is achieved by keeping the AAVS 0.5 beam 
pointed in a fixed direction while tracking the calibrator source 
with the MWA. Fig. [5] shows the X-polarization beam for 
such a measurement, where the MWA tracks Hydra A over 
a 4 hour period and the AAVS 0.5 beam remains pointed at 
Az=0°, El=75.4°. For each snapshot and channel, calibration 
is performed and A/T is calculated for the AAVS 0.5 beam at 
the direction of Hydra A. The data is normalized to the beam 
maximum. For clarity, only the 2% moisture case is shown, 
as the difference with the 10% case is not significant. 

Fig. [5] is a slice through the AAVS 0.5 beam pattern corre¬ 
sponding to the Az, El trajectory of Hydra A; the simulated 

4 http ://www. feko. info/ 
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Fig. 4. Measured A/T for a 2 minute observation of Hydra A starting at 
22-May-2014 17:30:32 at the MRO, X-polarization (East-West arm) top, Y- 
polarization (North-South arm) bottom. Each data point is a 40 kHz channel 
measurement. The AAVS 0.5 array pointing is Az=0°, El=75.4°. 


results in Fig. [5] are for the same trajectory. The solid orange 
curve simulates the real-world scenario where the sky “seen” 
by the AAVS 0.5 beam varies as a function of time, inducing 
changes in T sys . The dashed curve is simply the beam pattern 
(i.e. no sky is present). 

Towards the end of the observation (~23:10), the sidelobe 
level from the model sky curve is decreasing relative to the 
“no sky” curve; this is caused by a reduction in A/T due to 
increasing sky noise. However, the small difference between 
the two beam simulations is encouraging, indicating that the 
measured beam pattern (which intrinsically varies with time) 
is representative of the actual beam pattern, at least for the 
sky near Hydra A. 

We also estimate errors in the beam pattern introduced due 
to the analog RF beamforming system. Our estimate of the 
errors is calculated from measured phase (cr^ = 0.069) and 
amplitude (cr a =0.10) errors. The calculation for the ensemble 
mean beam pattern (sidelobe level) that incorporates these 
errors is well-known ED-ED; for our results there is no 
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Fig. 5. The AAVS 0.5 X-polarization beam pattern at 220 MHz for Az=0°, 
El=75.4° pointing. Hydra A was continuously observed on 21-Mar-2014 with 
2 minute snapshots and 5.12 MHz bandwidth. Each data point is a 40kHz 
channel measurement. The black curve is the same data, frequency-averaged 
for each snapshot. The oranges curves are simulated (“error-free”) beam 
patterns for 2% moisture in the presence of model sky (solid) or no sky 
(dashed). The purple curve is the 2(7 uncertainty on the beam pattern. 


significant difference to the error-free beam pattern. We plot in 
purple the estimated ±2a error on the mean pattern, where a 
is approximated using the formula in fl7| for large N. These 
curves show that main beam and most of the sidelobes are 
within 2cr of the mean beam pattern; the exception being the 
first sidelobe at 20:38, which is ^ 4a from the mean beam 
pattern. 

IV. Conclusion 

This work demonstrates far-field sensitivity and beam pat¬ 
tern characterization of a low-frequency aperture array. The 
measurements are in-situ, both in terms of physical location 
and intended astronomical use. The interferometric observa¬ 
tion of astronomical sources using the MWA radio telescope 
and the AAVS 0.5 array under test results in measurements cal¬ 
ibrated to an absolute scale. The sensitivity measurements and 
their frequency-dependent trends are consistent with simulated 
results. The measured and simulated beam patterns generally 
show good agreement as to the location and height of the 
sidelobes within the error. This type of characterization will 
be essential for verifying the performance of future prototype 
low-frequency aperture arrays deployed at the MRO. 
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